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SUMMARY 

i  .  1 

We  derive  a  formula  for  the  expected  number  of  crossings  of  an 

A 

arbitrary,  possibly  discontinuous,  curve  in  a  time  interval  (0,T) 
a  continuous  non-stationary  normal  process.  It  is  chown  that  a 
formula  analogous  to  that  of  M.  R.  Leadbetter  and  J.  D.  Cryer  [3] 
holds  under  conditions  more  general  than  they  considered. 


by 


Introduction 


We  consider  the  expected  number  of  (non-tangential)  crossings  of  an 
arbitrary,  possibly  discontinuous,  curve  in  a  time  T  by  a  continuous 
non-stationary  normal  process.  It  is  shown  that  a  formula  analogous  to  that 
derived  by  M.  R.  Leadbetter  and  J.  D.  Cryer,  holds  under  conditions  more 
general  than  they  considered  in  f 3 ] .  The  method  used  is  that  developed  by 
N.  Donald  Ylvisaker  in  [4],  For  the  sake  of  simplicity  we  assume  that  the 
process  and  the  curve  have  been  simultaneously  normalized  so  that  wc  may 
consider  crossings  of  a  curve  v(t)  by  a  continuous  process  with  mean 
function  0  and  with  variance  1.  It  is  shown  that  if  v(»)  is  of 
unbounded  variation,  then  the  expected  number  of  crossings  is 

The  proofs  are  somewhat  more  intuitive  than  those  of  [3]  in  that 
they  proceed  directly  from  a  calculation  of  the  probability  involved 
instead  of  an  interplay  of  Dirac  delta  functions  and  the  process 
derivative. 

The  formula  for  the  expected  number  of  crossings  is  a  sum  of  two 
parts;  the  first  is  just  an  integral  of  Rice's  formula,  the  derivation  of 
which  is  independent  of  the  second  part  and  follows  Ylvisaker,  14 J.  For 
v(*)  =  constant  this  gives  a  slightly  more  general  form  of  Ylvisaker's 
proof  of  the  necessary  and  sufficient  conditions  for  Rice's  formula. 

The  second  part  of  the  sum  represents  the  expected  number  of  crossings 
of  the  vertical  component  of  the  curve.  If  the  curve  is  a  step  function. 
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the  second  part  gives  the  sum  of  the  probabilities  of  crossing  the  vertical 
portions  of  the  curve.  Examination  of  this  part  can  give  information 
about  the  error  involved  if  we  replace  the  curve  with  a  step  function, 
(thereby  simplifying  the  necessary  calculations). 

As  a  corollary  we  get  a  simplified  form  of  the  result  of  Leadbetter 
and  Cryer  under  conditions  somewhat  more  general  than  they  consider.  The 
simplification  in  form  is  essentially  a  result  of  the  prior  normalization 
of  the  process  and  curve  to  get  a  process  with  variance  function  1.  If 
this  normalization  is  made,  the  hypothesis  in  [3]  that  fl<>(t,t')  is  con¬ 
tinuous  at  diagonal  points  is  equivalent  to  our  hypothesis  that  the  o(t) 
of  Theorem  1  is  continuous. 

The  expected  number  of  crossings  of  a  curve  of  unbounded  variation  has 
little  practical  application,  but  the  fact  that  it  is  always  +">  for 
continuous  Gaussian  processes  does  give  information  about  the  behavior 
of  such  process;  and  for  non-Markovian  processes  very  little  such 
information  is  available  in  the  literature. 


Notation 


Given  a  point  function  v(*)  of  bounded  variation,  let  v i • 1  be 
the  measure  determed  by: 

v{(o,t]}  =  .  (t)  -  (o) . 

Let  the  measure  j  v ) { • }  be  the  total  variation  of  Let 

|  v | ( • )  be  the  point  function  determined  by  j  e ' { • } : 

|  V  |  (t)  *  I v I {<o,t]}. 

We  shall  say  that  X(t)  has  ctossed  the  curve  (t)  in  the  time 
interval  T  if  there  exists  tj,  t?  in  T  such  that 
(X( t  j )  -  v(tj))  •  (X(t2)  -  v(t?))  <  0.  The  number  of  crossings  of  the 

curve  is  then  the  maximum  number  of  disjoint  time  intervals  in  which 
X(*)  crosses  v ( • )  • 

Theorem  1; 

Suppose  X(t)  is  a  non-stationary  Gaussian  process  with  mean 
function  0  and  variance  function  1. 

We  assume  that  the  correlation  function  p(tj,t>)  =  EX(ti)X(t2) 
has  the  spectral  form 

P ( 1 1 , t2 )  =  f  cos  2n\(t^-t])F(tltd>) 

—00 

with  the  extended  real  valued  continuous  on  (0,T)  second  spectral 

00 

o2(t)  =  J  (2irX)2F(t,d'). 


moment 


Let  N(T)  denote  the  number  of  crossings  of  a  fixed  curve  \  (t) 
by  the  process  X(t). 

Then : 


The  first  integral  in  (1)  is  just  the  integral  of  Rice's  formula 
for  the  number  of  crossings  of  a  fixed  level;  it  does  not  involve  /(•). 


If  v(*)  is  constant,  then  the  integral  with  respect  to  j  v  j ( • }  vanishes. 

On  the  other  hand,  if  '>(•)  is  a  step  function,  then  the  second 
integral  is  the  sum  of  the  probability  that  X(t)  crosses  the  various 
vertical  components  of  v(*).  In  this  sense  it  may  be  said  that  the 
first  term  represents  the  expected  number  of  crossings  of  the  horizontal 
component  of  v(*)  and  the  second  term  represents  the  expected  number 
of  crossings  of  the  vertical  component. 

Proof  of  Theorem  1: 

Let  d  *  { t ^ :  0  <  tj  <  t2  <  <  tn^  <  1}  be  a  collection  of 

points  of  the  interval  (0,T).  Let  y^(t)  be  the  stochastic  process 
whose  graph  is  a  series  of  straight  lines  with  vertices  at  (xft^)  -  '■’ ( t ^ ) ,  tj). 
Let  N^(T)  be  the  number  of  axis  crossings  by  y^(t). 

We  shall  be  interested  in  taking  limits  over  some  fixed  sequence  of 
collections  Id^fj  having  the  property  that  d^^  contains  d^,  and 
lim  max  (t^+^  -  t^)  *  0. 

*■“  VWdi 

For  the  portion  of  the  proof  dealing  with  v(*)  of  bounded  variation, 

oo 

we  further  assume  that  U  d,  contains  all  the  discontinuity  points  of 

1  i 

v(*) . 

Let  N,(T)  be  the  number  of  axis  crossings  by  y.(t),  0  i  t  i  T. 

d  a 

Then : 

N(T)  >  N  ,(T) 
a 


and  if  v(*)  is  of  bounded  variation 
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lim  N  ,(T)  -  N(T)  a.s. 

d 

hence 

lim  EN,(T)  -  EN(T) . 
d 

Taking  expectations: 

EN^CT)  *  Pf{yd(t)  has  an  axis  crossing  in  ( t ^ , t ^ 

-  Zpr<X(ti+i)  >  v(ti+1)}  -  Pr{x(t1+1)  >  v(ti+1).  Htt)  >  v(ti)} 

+  Pr{x(ti)  >  V(ti)j  -  Prix(ti)  >  v(ti),  X(ti+1)  >  v(ti+1)i- 

Let  L(h,k,P)  be  the  probability  that  Z  >  h,  W  >  k,  when  Z 
and  W  have  a  bivariate  normal  distribution  with  means  0,  variances  1, 
and  correlation  P. 


Then 


ENd(T)  -  £|L(v(t1+1),  v(ti+1),  1)  -  L(v(ti+1),  v(ti>, 

+  L(u{tp  i  v(tj)  (  1)  “  L(v(t^)  ,  •  P  ( )  ) 

Z|w?L(V(ti+l)*  "  L(v(ti+1)’  v^ti+l^’  p^ti+l*ti^^J 

(ti+l  '  ti) 


+  Z 


+  xi 


It ,,,-t  [L( v( ti) ,  v(ti),  1)  -  L(v(ti),  v(ti),  P(ti+1»t1))l[  *  (ti+1  -  tj) 


i+1  i 


r  V(W 


|v|(ti+1)-M(tl) 


/  / 

Lv(t1)  v(ti) 


^(x,y,P(tj+1,t1)dxdy 


*  (M(ti+1)  -  |v|(t4)) 


III 

J  nd1J(t)dt  +  /  n*2)(t)dt  +  /n*3*  (t)  |v|  (dt) 


where  n^Ct),  n^(t),  n^U)  are  the  step  functions  represented  in  the 
braces  in  the  above  sums. 


The  proof  of  the  theorem  follows  from  lemmas  1,  2,  3,  and  4. 


To  simplify  the  notation  let  t  be  an  element  of  the  collection 


(and  hence  an  element  of  all  subsequent  d^) ,  and  let  :  be  the  de¬ 
creasing  difference  from  t  to  its  successor  in  d^.  We  will  abbreviate 
p(t,t+r)  by  p(t). 

Lemma  1.  Under  the  conditions  of  Theorem  1,  if  v(»)  is  of  bounded 

variation,  then: 

T  T 

lim  /  nJ1J(t)dt  *>  ~  f  e_!sv2(t:)[7  (2tt>)  ?  F  ( t ,  d  X  )J  *^dt 

0  "  0  0 

in  the  sense  that  if  one  side  of  the  equality  is  then  both  sides 

are. 


-8- 


Letting 


-y- 


we  have 


.  [arC  <T>]  =  a'5(t) 


1  3  t"A  , 

+  2  *  4  *  - 5 - N  0 . 


Therefore 


-  Cju  _ 

>-r.c..cosJ,(T)j  <  2J(l=j 2.cos2,  ■:|F(t>d.) 

*—  — 


-  X 

-  —  1  1m  *1 

,  f  /  1-C0s2'1>.  T  \ 

=  lim  i 

J  [  ,  ) 

J  TiO 

•'o'  ' 

in  the  sense  that  if  one  side  of  the  equality  is  4^  then  both  sides 
are  4°°. 

1-cos2tt>,t 

But  the  non-negative  function  - -  is  bounded  above  by  and  con- 

•t 

T  ‘ 

verges  to  2tt2A?.  It  follows  that 


K(i)  Irz  I ( 2 ^ A ) •’ F( t , d o  ' 


oc 

C(t)  =  e”'2J  (t  )  /(2nA)-F(  t,d\) 


lim  K( 
TlO 


in  the  sense  that  if  one  side  of  these  equalities  is  +“*  then  both  sides 
are  +°°. 

The  statement  of  Lemma  1  follows  since  v(t)  =  v(t+)  a.e. 

I 

The  Proof  of  Lemma  2  follows  in  the  same  way. 


Proof  of  Lemma  3.  If  t^  is  a  point  of  discontinuity  for  v(*), 

it  is  a  point  of  positive  mass  for  |  v  1  { • }  and  we  must  do  better  than 

( 3) 

to  find  the  a.e.  limit  of  n,  (t)  at  such  points.  For  the  sake  of 

d 
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convenience  we  will  assume  that  t  =  t^  is  a  point  of  left  continuity  of 
v(*);  the  necessary  modifications  are  recognizable  if  this  is  not  the 
case. 


Letting  t  =  t1+J  -  t^  p(t)  =  p(t,t+x): 


(3)m  =  1 

d  Tv  (t+x)  -  v  (t) 


v(t+i)  v(t+x) 

/  /  •<*  ,v ,p ( i ) )  dxdv 

v(t)  v ( L) 


1 

]v~](t+x)  -  I  V  I  (  t ) 


v(t+x) 

/ 


v(t) 


v(t+;)-p ( i  )x 

Vi-p-'(i) 

$(x)dx  I  i£(y)dy. 


v(t)-n(i)x 
v/l-p’(  0 


/  o\ 

(Tliat  n^  ^(t)  is  bounded  above  by  a  |v|{»}  integrable  function 
follows  if  we  take  1  as  an  upper  bound  for  the  integral  over  y  and 
<P(o)  as  an  upper  bound  for  the  integrand  of  x.) 


Substituting 


x  =  pv(t+x)  +  (l-p)v(t),  0  <  p  v  1, 


let 


1  t) 

M(t)  =  |  vj’[t+I)  I  Jvflt)  J  +  ( 1-p)  v( t))dp  •  ^(v)dy 

a 

1  avb 

=  -  jv|(t)  ■'  i( Pv( t+T )  +  (l-p)v(t))dp  •  f  y)dy 


aAb 


where 
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a  =  a(o.T)  =  v.(.t)  ~  LL-liEHjt±ll  +  lizniiilU 

/l  -  P‘‘(t) 

b  =  b(p>T)  =  ^(1+^)  Lp- (t+T)  ±JLj~£iLLLll 

/l  -  p2(x) 


To  evaluate  lim  a(p,x),  write 
T  10 


V 1  -  o2(0 


^v(t  +  x)-v(t)  j  . 


+  o(  l)  . 


A  -  C-'(t) 

oo 

Since  0  <  f  (27tA)^F(t,d\)  <  •+“>,  t  (•)  is  twice  differentiable,  and 


we  may  evaluate  lim  t />/ 1  -  p2(t)  by  squaring  and  applying  L’llospital ' s 
rule : 


lim 


1  -  P2(i)  -p"(o)  o-'(t) 


therefore 


lim  a(p,x)  =  -pv  (t)/o(t)  uniformly  in  p,  0  1  p  i  1, 

T*0 


Similarly 


lim  b(p,T)  =  (l-p)v  ( t) / o ( t) 
TfO 


uniformly  in  p,  0  i  p  :  1, 


Thus 


avb 


(l-p)|v+(t)|/o(t) 

/  ^(yjdy  =  f  $(y) dy . 

J '  J  i  +  . 


aAb 


-P'1  v  (t)  |  /a  ( t ) 


*• 


It  follows  from  the  orthogonal  decomposition  of  a  signed  measure 


that 


lim 
t  10 


N(t+0  -  y(t)  j 

vj  (  t  +  T  )  -  I  V  [  (  t  ) 


1. 


Hence 


M(o)  = 


1  (l-P)|v+(t)|Mt) 

J  4>(pv(t+)  +  ( 1-p) v ( t) ) dp  I  c(y)dy 


-p|v+(t)  |/-(t) 


1  p I v  (t)  j h  (t) 

J  $(pv(t  )  +  ( 1-p) v  ( t) )  dp  /  ,"(y)dv 

r\ 


~P  I  v  (t)j/o(t) 


I 


If  v(»)  Is  of  Unbounded  Variation 

Lemma  4.  Suppose  v( • )  is  a  function  of  unbounded  variation  and 
o(t)  is  bounded  or  is  extended  real  valued  continuous. 

Then 


1  im 

max(  t^+^-t . )  *0 


Ul 


v(ti+l> 


(ti> 


:(x,y,.  (t.+1,t .  ))dxdy  =  +■> 


Corollary .  It  follows  from  Fatou's  lemma  and  the  inequality  N('l)  _ 
that  if  v(*)  is  of  unbounded  variation  and  o(t)  is  continuous  or  is 
bounded  on  (0,T),  then  EN(T)  =  +nc>. 

Proof  of  Lemma  4.  As  in  Lemma  3,  we  have: 


v(ti+l)  .•(ti+l) 


M(d)  =  l  f  f  : (x,y ,c (ti+]  ,t . ))dxdv 


v ( 1 1 )  v(t.) 


1  b 

X  (  v  < 1  i-n>  "  /  :(Py(ti+j)  +  ( 1-p)  o  ( t  ) )  dp  fi(x)dx 

0 


a 
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where  a,b  are  as  in  Lemma  3.  We  may  take  a  more  conservative  range 
of  integration  by  setting 


p|v(t  )-v(t  )  1-;  (t  t  ) 

^ ' uv -'■‘w + 


"i+l  i' 


then 


c.  (p) 


M(d)  >  2^|v(ti+1)  -  v(ti)|/ idpj1  <dx. 


Let  0  =  max  |v(t)|,  0  <  t  <  T,  and  fix  t,  0  «•  r  «■  1;  then 

c .  (p) (1-e) 

M(d)  >  c  $(v)£  |  v(ti+1)  “  v ( tA)  |  •  ^  ,"(x)dx. 


Define  the  set  Sj  depending  on  d  =  {  t ^ /  as  follows: 


!v(tH-l>  -  "(ti>: 

/I  -  o(tJ+1.tt) 


7fc•)(5+;,, 


=  { i :  ci(l-0  2  5}. 


Then  for  i  e  sd  we  liave 

c^p)  (1-c)  b 

J  (t(x)dx  1  J  |(x)dx  =  6, 

'o  0 

M(d)  j.  c  •  6  •  4 ( v) [  v  ( t  !  • 

Sd 

If  v(*)  is  a  function  of  unbounded  variation,  then  this  sum  over 
all  i  is  unbounded  for  any  sequence  of  d  having  the  property  that 


lim  max  (t^^-t^)  1 

t^cd 


thus  Lemma  4  follows  from  Lemmas  5  and  6. 


I 


rnMKmmmtm 


* 
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Le mma  5.  The  sequence  =  <t^,t7,...[  may  he  picked  so  that 


}_'  y/l  -  o'  (tj+j»tj)  is  bounded  above  for  all  i;  and  max  (t^+^-t.)  •  0. 


Proof , 


n/T~-  px(tH1,t.) 

Cj  +  1  "  ‘j 


(l-cos2- ' (t ,-t  .)) 

(l+c(t  t  .))[ - L~ - L—  F  <dO 

J  A»  J  J  It  -t  V  cj 


H. 


“j+l'V 


F  Id- 


11 


If  o(»)  Is  bounded  on  (o,T)  we  are  through.  If  it  is  continuous 
and  integrable,  then  for  a  given  i,  let  the  t.  j  =  in  d. 

be  picked  so  that  o ( t j )  =  min  o(t),  t  in  the  interval 
( J / 2 1 ) T  I  t  :  (j+l/2l)T.  Then 

T 

}'  /l  -  (tj+1,t  ^  <  .Ht  j+1-tj)  :  /  (t)dt . 


I 


Lemma  6.  Let  ajn^  40  be  a  sequence  such  that 


L  Va!n)  = 


n  —  l 


Let  b|n^  >  0  be  a  sequence  such  that 

r 

^  bfn^  K,  some  fixed  K. 

—  l 

I 

For  arbitrary  fixed  M,  let  S  =  ^i:  ain^/bfn^  '  m}  . 


Then : 


L  l  a!n)  - 
nLS  i 
n 
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Proof  . 
bounded: 


wnmm.t  Aim.miiGmirxsgsxmittiZtSCSe  S)  mMKt  _ 

I, 


Let  Sc  =  li:  a. /b . 


We  show 


V  a(n) 

s‘  1 


is  uniformlv 


)  a,(n) 
:c  i 

3 

n 


n 


n 


I 


Corollary.  (Leadbetter  and  Cryer)  Suppose  that  the  function  v(*) 

(and  hence  the  function  |  v |  ( • )  and  the  measure  •  v  |  { * } )  are  absolutely 
continuous  with  respect  to  Lebesgue  measure.  Then  the  formula  (1)  of 
Theorem  1  reduces: 


T 

r  ,  v*(t)/a(t)  , 

(2) 

EN(T)  =  J 

e  ^  ^  +  o(t)f(v(t)) J  (2<f(y)-l)dy 

0 

0 

Proof.  In  lieu  of  the  absolute  continuity  of  j  v | { • }  we  may  rewrite 
(1)  as  an  integral  with  respect  to  Lebesgue  measure;  consider  the  integrand 
of  the  second  portion  of  that  integral: 

n(t)  *  c(o,v,t)  •  J  v * ( t ) | 

1  p  !  V  *  ( t )  j  /o  (  t ) 

*  2  j  v * ( t ) | $ ( v ( t ) ) f dp j  <(y)dy  , 

0  0 


1  p(v' (t)/o(t)) 

n(t)  =  2o  ( t)  <f  (v(  t)  )  ( v '  ( t) /o  ( t) )  J  dp/  ■f(y)  dy. 

0  0 


l 


Writing  £  =  v'(t)/o(t),  and  considering  n(t,£)  =  n(t): 

n(t,0)  =  0 


and 
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— -  (t)(24(;)-i). 

Therefore, 

n(t,0  -  (2i(y)-l)dy. 

0 

I 


Hie  aovollavj  is  unlike  the  result  of  [31  only  in  that  we  have 
assumed  v(*)  absolutely  continuous  instead  of  continuously  dif fcrentiabl 
i f  we  consider  a  process  with  variance  function  1,  then  the  hypothesis 
of  [3]  that  F  12(t , t * )  exist  and  is  continuous  at  diagonal  points  is 
exactly  our  requirement  that  o(t)  exist  and  is  continuous.  To  show  that 
the  formula  of  [3]  reduces  to  (2)  when  the  variance  function  is  1,  note 
the  typographical  error  on  p.  510: 

it  follows  (adding  a  circumflex  to  denote  the  notation  of  (3J) 
from  Lemma  3  of  [3]  that 

n(t)  =  (m'-ypni)  / [ y  ( l-v‘ )  Jj. 

If  the  variance  is  unity,  the  following  simplifications  may  be  made: 

0"  x  1 
YJ  -  aJ(t) 
c:  «  0. 

With  these  substitutions  the  formula  (2)  of  [3]  may  be  reduced  to  (2) 
above  by  the  same  method  of  dif ferentiating  with  respect  to  m'(t)/o(t). 
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